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The search for di-Higgs production at the LHC in order to set limits on Higgs trilinear coupling 
and constraints on new physics is one of the main motivations for the LHC high luminosity phase. 

Recent experimental analyses suggest that such analyses will only be successful if information from 
a range of channels is included. We therefore investigate di-Higgs production in association with two 
hadronic jets and give a detailed discussion of both the gluon- and weak boson fusion contributions, 
with a particular emphasis on the phenomenology with modified Higgs trilinear and quartic gauge 
couplings. We perform a detailed investigation of the full hadronic final state and find that hhjj 
production should add sensitivity to a di-Higgs search combination at the HL-LHC with 3 ab _1 . 

Since the WBF and GF contributions are sensitive to different sources of physics beyond the Stan¬ 
dard Model, we devise search strategies to disentangle and isolate these production modes. While 
gluon fusion remains non-negligible in WBF-type selections, sizeable new physics contributions to 
the latter can still be constrained. As an example of the latter point we investigate the sensitivity 
that can be obtained for a measurement of the quartic Higgs-gauge boson couplings. 


I. INTRODUCTION 

After the Higgs boson discovery in 2012 0 and subse¬ 
quent analyses of its properties [4], evidence for physics 
beyond the Standard Model (BSM) remains elusive. Al¬ 
though consistency with SM Higgs properties is expected 
in many BSM scenarios, current measurements do not 
fully constrain the Higgs sector. One coupling which is 
currently unconstrained and has recently been subject of 
much interest is the Higgs self-interaction ~ r /, which is 
responsible for the spontaneous breaking of electroweak 
gauge symmetry in the SM via the potential 

V(HiH)=n 2 H*H + ri(HiH) 2 , ( 1 ) 

with fi 2 < 0, where H = (0,u + h) T /\pi in unitary 
gauge. The Higgs self-coupling manifests itself primar¬ 
ily in a destructive interference in gluon fusion-induced 
di-Higgs production M through feeding into the trilin¬ 
ear Higgs interaction with strength Asm = trih\J J?/2 = 
gm 2 J(Aniw) in the SM. The latter relation can be al¬ 
tered in BSM scenarios, e.g. the SM coupling pattern 
can be distorted by the presence of a dimension six op¬ 
erator ~ (WH) 3 , and di-Higgs production is the only 
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channel with direct sensitivity to this interaction [6j. A 
modification solely of the Higgs trilinear coupling, which 
is typically invoked in di-Higgs feasibility studies, is pre¬ 
dicted in models of [i 2 -\ess electroweak symmetry break- 
ing, e.g. 0. 

After the Higgs discovery, analyses of the di-Higgs final 
state at the high-luminosity LHC and beyond have ex¬ 
perienced a renaissance, and di-Higgs final states such as 

the bb 77 li-K , bbr + T~ [lTI4T3j| . bbW + W~ [TTI.[l2lITij 

and bbbb |ll], EH channels have been studied phe¬ 
nomenologically, often relying on boosted jet substruc¬ 
ture techniques [Iff (see also an investigation na of rare 
decay channels relevant for a 100 TeV collider). Recent 
analyses by ATLAS and CMS [lU, [l9[ have highlighted 
the complexity of these analyses and the necessity to ex¬ 
plore different production mechanisms to formulate con¬ 
straints on the Higgs self-interactions in the future. This 
program has already been initiated by feasibility anal¬ 
yses of the hhj , hhjj and tthh production modes in 

Refs. E3,[HlMll. ' 

Di-Higgs production in association with two jets is a 
particularly important channel in this regard since this 
final state receives contributions from the weak boson 
fusion (WBF) production mode. The phenomenological 
appeal of the WBF mode is twofold. Firstly, the weak 
boson fusion component of pp -A hhjj is sensitive to 
modifications of the gauge-Higgs sector [2(1 EH, which 
can lead to large cross-section enhancements. Secondly, 
the QCD uncertainties for the WBF topologies are known 
and under theoretical control [24| [25| , such that a search 
for BSM electroweak-induced deviations is not hampered 
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FIG. 1: Maximum Higgs and jet transverse momenta for gluon fusion-induced hhjj production, including the ratio of full 
theory to the effective theory calculation for three different values of the Higgs trilinear coupling A. 



FIG. 2: Sample Feynman diagrams contributing to pp —> 
hhjj via gluon fusion. 


by QCD systematics. This situation is very different from 
QCD-induced production [2(|, and can be attributed to 
the particular phenomenology of WBF-like processes [23- 

m 

However, an additional source of uncertainty that was 
neglected until recently [2C| is the correct inclusion of 
the gluon fusion contribution to pp —> hhjj analyses. 
In contrast to single Higgs phenomenology, the correct 
inclusion of massive fermion thresholds is crucial to a 
reliable prediction of QCD-induced pp —I hhjj fl^ . 

Given that the cross sections in WBF hhjj production 
are very suppressed compared to WBF hjj production 
(the WBF hhjj cross section is ~ 750 times smaller), we 
have to rely on the dominant hadronic Higgs decay modes 
to be able to observe this final state. This rules out one of 
the most crucial single Higgs WBF selection tools - the 


central jet veto [29[. The observation of WBF-induced 
pp -A hhjj production is further hampered by the top 
threshold in the QCD-mediated process. Since the top 
threshold sets the scale of the di-Higgs subsystem, an 
analysis that tries to retain as many low pp Higgs bosons 
as possible leads to a QCD contribution that dominates 
over the WBF component when minimal WBF-like cut 
requirements are imposed [2C|. 

In this paper we extend the discussion of Ref. [2(J in 
a number of directions. We first perform a detailed com¬ 
parison of EFT-approaches to QCD-mediated pp —> hhjj 
against a calculation keeping the full mass dependencies 
of top and bottom quarks in Sec. El We compare the 
QCD-induced pp -A hhjj phenomenology to the WBF 
signature in Sec. IIIII before we discuss general approaches 
to isolate the signal from the dominant top backgrounds 
in a hadron level analysis in Sec. IIVI This sets the stage 
for a discussion about the prospects to isolate the WBF 
and GF components in Secs. IIV Al and HVBl followed by 
a study on constraining VVhh couplings using the WBF 
induced signal in Section IIV Cl We focus on collisions 
with 14 TeV throughout. 


II. THE GLUON FUSION CONTRIBUTION 
A. Finite top mass effects 

It is well known that effective field theory approxima¬ 
tions in the mt —> 00 limit cannot be invoked to study 
di-Higgs final states at colliders in a reliable way due to 
the effects of top-quark threshold H, H3- Further, the 
breakdown of the m* —> 00 approximation is worsened 
in the presence of additional jet emission [l2|, ETJ . Finite 
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A 

0 • Asm [fb] 

1 ■ Asm [fb] 

2 • Asm [fb] 

GF 

10.73 

5.502 

2.669 

WBF 

4.141 

2.010 

0.9648 


TABLE I: Cross section normalisations for the GF and WBF 
samples at 14 TeV, for details see text. The WBF normalisa¬ 
tion follows from [24j] and includes higher order QCD effects. 


TO* effects must therefore be considered for all QCD di- 
Higgs production channels, which will be required to set 
the best limits on the Higgs self-coupling or formulate a 
realistic estimate of the GF contribution in a WBF-like 
selection. 

The computational challenges in QCD-mediated hhjj 
production are significant, with the gluon-fusion channels 
particularly time consuming even when using state-of- 
the-art techniques. The standard method of simulating 
a differential cross section from unweighted events is not 
feasible in this case, and we instead use a reweighting 
technique that is exploited in higher order calculations 
and experimental analyses (see e.g. [32}). 

We generate GF hhjj events by implementing the 
relevant higher dimensional operators in the mt —> oo 
limit obtained by expanding the low-energy effective the¬ 
ory n 

C eS = -\^G; u G a ^log(l + h/v) (2) 

in MadEvent v 5.1 [34j using the FeynRules/Ufo 135 
framework.* This allows us to sample a weighted set of 
events that we subsequently feed into our analysis solely 
depending on their final state kinematics. If an event 
passes the selection requirements of a certain search re¬ 
gion, we correct for the full mass dependence using a 
reweighting library based on GoSam package Q at 
this stage. The reweighting employs exactly the same 
matrix elements used for the event generation and the 
trilinear coupling is steered through a modification of 
the GoSam matrix element, i.e. variations of the tri¬ 
linear coupling are part of the reweighting. A selection 
of Feynman diagrams which contribute to the gluon fu¬ 
sion signal are illustrated in Fig. [2j The GoSam code 
used for the reweighting is based on a Feynman diagram¬ 
matic approach using QGRAF §3 and FORM [3g for 
the diagram generation, and Spinney [3iJ, Haggies [4G] 
and FORM to write an optimised fortran output. The 
reduction of the one-loop amplitudes was done using 
Samurai |4lj| . which uses a d-dimensional inte gra nd level 
decomposition based on unitarity methods [421 ] . The 
remaining scalar integrals have been evaluated using 
OneLoop [43j . Alternative reduction techniques can be 


*The effective theory implementation can be modified in the sense 
that only one effective vertex insertion is allowed. This is gives 
only a mild ~ 10% effect in the tail of the distribution, and is not 
relevant for an order one EFT/full theory rescaling, see below. 


used employing Ninja J44| or Golem95 [45[. To vali¬ 
date the reweighting procedure we regenerated the code 
that has been used in j2d ] with the improvements that 
became available within GoSam 2.0, in particular im¬ 
provements in code optimisation and in the reduction 
of the amplitudes. For the reduction we used Ninja, 
which employs an improved reduction algorithm based 
on an Laurent expansion of the integrand. This leads 
to substantial improvements in both speed and numeri¬ 
cal stability compared to the previous version. We com¬ 
bined the code with a phase space integration provided 
by MadEvent [4f|. Further substantial speed-up has 
been obtained by Monte Carlo sampling over the helic- 
ities rather then performing the helicity sum. This en¬ 
abled us to perform a full phase space integration and 
we found full agreement within the statistical uncertain¬ 
ties between the result obtained from reweighting and 
the result from the full phase space integration. 


B. Phenomenology of QCD-mediated hhjj 
production 


Top thresholds are particularly prominent in the di- 
Higgs invariant mass distribution, which is thus well 
suited to benchmark the relation of the finite m t limit 
to the effective theory of Eq. ©• Other observables con¬ 
structed from the six particle final state are also relevant 
when performing a targeted phenomenological analysis, 
and we discuss both these and the phase space dependent 
parton-level reweighting in detail in the following. 

In Figures [Q [Sj and [I] we show a selection of hhjj fi¬ 
nal state observables for inclusive cuts pr,j > 20 GeV 
and |? 7 j| < 5, no cuts on Higgs bosons are imposed. We 
label Higgs bosons and jets according to their hardness, 
i.e. PT,hi > PT,h2 and Ptji > Ptj 2 - The cross sections 
are given in Tab. Q] The inclusive gluon fusion cross sec¬ 
tion is about 2.5 times larger than the WBF cross section 
approximately independent of the value of the Higgs tri¬ 
linear coupling. 

As previously established in SEH3 the di-Higgs 
system is badly modelled by the effective theory which 
under- and overshoots the full theory cross section at low 
and high momenta respectively. For pp —> hhjj this is a 
qualitatively similar behaviour compared the pp —> hh(j) 
production: The rrihh distribution is the crucial observ¬ 
able which parametrises the finite top quark mass effects. 
The EFT describes the low maximum transverse Higgs 


momenta Pt,m reasonably well, as shown in Fig. 1(a) 


The jet emission on the other hand integrates over a con¬ 
siderable range of rrihh , and the ratio of full theory vs ef¬ 
fective theory smaller than one for the finite mt limit pro¬ 
duces a smaller integrated cross section than the mt —> oo 
limit for the jet kinematics. 

Considering just the dijet system in Fig. [3l we observe 
that the jet kinematics is not severely impacted by the 
reweighting procedure upon marginalising over the di- 
Higgs kinematics. The phase space dependence of the 
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FIG. 3: Invariant mass and pseudo-rapidity distributions of the jet system in QCD-mediated hhjj production. We show the 
effective theory and full theory results for three values of the trilinear Higgs coupling, applying only generator-level cuts of 
Pr,j > 20 GeV and \rij\ < 5. 
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FIG. 4: Invariant mass and pseudo-rapidity distributions of the leading jet and di-Higgs system in QCD-mediated hhjj 
production. We show the effective theory and full theory results for three values of the trilinear Higgs coupling, applying only 
generator-level cuts of pr.j > 20 GeV and |iy \ < 5. 


dijet invariant mass Fig. 3(a) is relatively mild aside from ported by the larger impact of the reweighting of mhh 


the total rescaling of the inclusive cross sections, and the 
ratio for the pseudo-rapidity distribution of the jets is 
nearly flat, Fig. 3(b) This is also true for the azimuthal 
angle difference A (f>jj. The angular distributions of the 
leading members of the jet-Higgs system ar e rela tively 
mildly impacted by the reweighting too Fig. |4(b)| This 
agrees with the mhh being the observable most sensitive 
to the top threshold (as in pp hh(j)), and is also sup- 


in Fig. 4(a) A reweighting based on mhh to correct for 


finite top mass effects suggests itself for future analyses 
as a time-saving approach with reasonable accuracy. 
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FIG. 5: Sample Feynman diagrams contributing to pp —> hhjj 
in weak boson fusion. 


III. THE WEAK BOSON FUSION 
CONTRIBUTION 


The weak boson fusion contribution to pp —> hhjj 
has received considerable attention recently and pre¬ 
cise hi ghe r-order QCD corrections have been provided 
in [24], [25j, 0 . Due to the sensitivity of the WBF con¬ 
tribution to both the trilinear coupling and the quartic 
VVhh (V = W, Z , 7 ), as shown in the Feynman diagrams 
in Fig. 0 weak boson fusion to two Higgs bosons can, 
in principle, provide complementary information about 
BSM physics which remains uncaptured in pp —> hh(j) 
and pp —f tthh [23j . 

We generate WBF samples with varying A using 
MadEvent v 4 [471 and normalise the cross section to 
NLO accuracy [24j. The WBF hhjj contribution shares 
the QCD properties of WBF hjj production [2jJ which 
means it shares the distinctive Arj(jl, j2) distribution 
shown in Fig. | 6 (a)| To produce the heavy di-Higgs pair 
we probe the initial state partons at large momentum 
fractions. This together with a colour-neutral f-channel 
exchange of the electroweak bosons [HI (see also 0) 
leads to energetic back-to-back jet configurations at large 
rapidity separation and moderate transverse momenta 
with a centrally produced Higgs pair. The production of 
an additional Higgs boson in comparison to single Higgs 
production via WBF leads to a cross section reduction by 
three orders of magnitude (see Tab. QJ in the SM. Such 
a small inclusive production cross section highlights the 
necessity of considering dominant Higgs decay channels 
such as h — > bb and h — > t + t~ and the non-availability 
of central jet vetos 0] as a means to control the back¬ 
ground and GF contribution in a targeted analysis as a 
consequence. 

The gluon fusion contribution is bigger by a factor of 
2.5 than the WBF component of hhjj production, how¬ 
ever, with increasing invariant di-Higgs mass the WBF 
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FIG. 6: The Ar/(jl,j2) distribution of the weak boson fu¬ 
sion contribution at parton-level (a) and the mhh distribu¬ 
tion of the weak boson fusion and gluon fusion contributions 
compared with correct cross section normalisation (b), both 
satisfying generator-level cuts of pr,j > 20 GeV and | r/j\ < 5. 


contribution is enhanced relative to GF production as a 
consequence of the suppression above the 2mt threshold, 
as shown in Fig. | 6 (b)[ 

Since we cannot rely on vetoing hadronic activity in the 
central part of the detector, a potential discrimination of 
GF from WBF needs to be built on the following strategy, 
which we will investigate in Sec. IIVI 

• To isolate the di-Higgs (WBF+GF) signal we can 
exploit the relative hardness of the di-Higgs pair 
which peaks around ~ 2 mt- Such hard events are 
less likely to be produced by (ir)reducible back¬ 
grounds. 
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Cut setup 

Base Selection [fb] 

GF Selection [fb] WBF Selection [fb] 

Normalisation* [fb] 

GF (A = 1 • Asm) 

1.396 x 10~ 2 

5.722 x 10 -3 

5.378 x 10“ 4 

4.013 x 10 1 

GF (A = 0 ■ Asm) 

2.562 x 1CT 2 

8.122 x 10 -3 

8.767 x 10“ 4 

7.831 x 10 _1 

GF (A = 2 ■ Asm) 

7.167 x 1CT 3 

3.906 x 10“ 3 

3.034 x 10" 4 

1.947 x 10 _1 

WBF (A = 1 ■ Asm) 

3.292 x 10“ 3 

4.999 x 10“ 4 

1.485 x 10“ 3 

1.466 x 10 _1 

WBF (A = 0 • Asm) 

7.706 x 10“ 3 

7.154 x 10“ 4 

2.820 x 10“ 3 

3.020 x 10 _1 

WBF (A = 2 ■ Asm) 

1.103 x 10” 3 

1.815 x 10“ 4 

3.912 x 10“ 4 

7.037 x 10“ 2 

ttjj 

5.712 

3.390 x 10“ 2 

1.801 x 10“ 2 

1.0130 x 10 4 

tth 

6.229 x 10~ 2 

7.047 x 10“ 3 

5.658 x 10“ 5 

3.862 x 10 1 

Zhjj 

5.118 x 10” 3 

1.278 x 10“ 3 

1.026 x 10“ 4 

4.737 x 10 1 

ZZjj 

1.171 x 10“ 3 

6.659 x 10" 5 

7.639 x 10“ 7 

2.257 x 10 2 

ZWWjj 

1.888 x 10” 5 

5.461 x 10“ 6 

2.039 x 10“ 7 

5.368 x 10" 1 

total background 

5.781 

4.230 x 10" 2 

1.817 x 10“ 2 

- 

S/B (A = 1 ■ Asm) 

1/335.1 

1/6.799 

1/8.983 


S/B GF f (A = 1 ■ Asm) 

1/414.3 

1/7.480 

1/36.55 


S/B WBF 1 * (A = 1 • Asm) 

1/1760 

1/96.06 

1/12.60 


S/V B (3 ab -1 , A = 1 • Asm) 

0.3930 

1.657 

0.8219 



* branchings included in normalisation 
1 considering only this as signal 


TABLE II: Cross sections for the two sources of signal, and backgrounds, after the various selections described in the text are 
applied, together with various measures of significance in the bottom four rows. 


• Focussing on large rrihh we can enhance WBF over 
GF by stringent cuts on the jet rapidity separation. 
This will also imply a significant decrease of QCD- 
dominated backgrounds. 

• By explicitly allowing central jet activity, we can 
exploit the colour correlation differences in WBF vs 
GF to further purify our selection. Since colour flow 
is tantamount to energy flow in the detector, event 
shapes are particularly well-suited observables for 
unravelling the colour correlations in the final state 
once the reconstructed di-Higgs pair has been re¬ 
moved^ . This strategy was first proposed for single 
Higgs analyses in 0 (see also |52jj. 

IV. TAMING THE BACKGROUND 

For our hadron-level analysis we shower our signal sam¬ 
ples with Herwig-I — b 0] and generate backgrounds as 
follows: ttjj , tth 1 Zhjj, and ZZjj with Sherpa 0], and 
ZWWjj with MadEvent v5. We find the dominant 
backgrounds to be ttjj and tth production^ for which 
next-to-leading order results are available 0, 0] and 
we use inclusive K factors K t ijj ~ 1 and K t ^ h ~ 1.5 to 
estimate the higher order contributions to these back¬ 
grounds. Higgs branching ratios are set to the val¬ 
ues agreed upon by the Higgs Cross Section Working 
Group 0. 


detailed discussion of event shapes at hadron colliders can be 
found in [50ll . 


We begin the hadron-level analysis implemented in 
Rivet 0 by recreating a base selections similar to [2(jjd 

1. ) We require two tau leptons using a two tau-trigger 

based on staggered transverse momentum selection 
cuts pr > 29, 20 GeV in \r/ T \ < 2.5 and assume a flat 
tau tagging efficiency of 70% with no fakes. 

Jets are constructed by clustering R = 0.4 anti -kr 
jets using FastJet 0 with prj > 25 GeV and 
M < 4-5. 

2. ) The two leading jets are 5-tagged with an acceptance 

of 70% and fake rate of 1% [0 in the central part 
of the detector | rjj\ < 2.5. We remove events if ei¬ 
ther of the two leading jets overlaps with a tau. Any 
additional jets which do not overlap with a tau are 
considered as potential “tagging jets”, of which we 
require at least two. 

3. ) As a final step of this base selection we require the b 

jet and tau pairs to reproduce the Higgs mass of 125 
GeV within ±15 and ±25 GeV respectively.§ 

The signal and background cross sections after these cuts 
are presented in the Base Selection column of Table HT1 
We find that the background contribution of ttjj domi¬ 
nates with tth also providing a larger-than-signal back¬ 
ground resulting in S/B ~ 1/300, making a study based 


tOur analysis has been validated with two independent implemen¬ 
tations. 

§A high mass resolution is a crucial cornerstone of any successful 
di-Higgs anal ysis to assure a minimum pollution of Z boson decay 
backgrounds 1131 . 
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FIG. 7: Shape comparison of Ar/(jl, j2) and rrihh distributions for our two sources of signal (GF and WBF), the dominant 
background tijj and the rest of the backgrounds (stacked scaled by relative cross sections), after the Base Selection of Section IWl 
has been applied. 


only on these selections extremely challenging. Since we 
only have ~ 40 expected gluon fusion and ~ 10 expected 
weak boson fusion events at 3 ab -1 luminosity, additional 
selections must also be careful to retain enough signal 
cross section to allow statistically meaningful statements 
to be made with a finite amount of data. 

Shape comparisons for the rapidity and dihiggs invari¬ 
ant mass distributions as motivated in the previous sec¬ 
tion are shown in Fig. [7] Indeed, as expected, cutting on 
the angular distance of the jets will serve to both purify 
towards a WBF-only selection at a reduced background 
rate. The dominant backgrounds are unlikely to produce 
a large invariant mass rrihh- However the WBF contri¬ 
bution, due to the lack of the 2m* threshold peaks at a 
considerably lower invariant mass, leading to significant 
decrease of the WBF contribution for a reasonably strong 
cut on rrihh, which is required to observe the hhjj signal 
at the given low signal yield, even at 3 ab -1 luminosity. 


A. Prospects to isolate gluon fusion 

We can extend the analysis outlined in Sec. IIVI with the 
aim to purify the selection towards the GF component.^ 
We make use of the hard Higgs candidates to greatly 
reduce the backgrounds by requiring nihh > 500 GeV and 
additionally require Ar)(jl, j2) < 5 to minimise the weak 
boson fusion contribution. The signal and background 
cross sections after these cuts are applied are presented 
in the ‘GF Selection’ column of Table [TT] 


^Following the analysis of [fi l i] . we can expect negligible interference 
between WBF and GF and which allows us to make this distinction. 


The total background is reduced by a factor of ~ 100 
while the gluon fusion contribution only is reduced by 
a factor of ~ 2.5 which allows for an encouraging 
S/V B ~ 1.7 with 3 ab -1 of data. The weak boson 

fusion contribution is also suppressed compared to GF 
which allows for a clean probe of the physics accessible 
in the gluon fusion contribution. 

B. Prospects to isolate weak boson fusion 

Similarly we can extend the analysis towards isolating 
the WBF component. Since it has slightly softer Higgs 
candidates we require rrihh > 400 GeV and Ar/(jl, ;/2) > 
5 to reduce both the gluon fusion and background contri¬ 
butions. The signal and background cross sections after 
these cuts are applied are presented in the ‘WBF Selec¬ 
tion’ column of Tabic (TT] 

The total background is reduced by a factor of ~ 300 
while three times more of the weak boson fusion contribu¬ 
tion is retained compared to the GF selection, resulting 
in S/\/~B ~ 0.8 with 3 ab -1 of data due to the large 
reduction in the gluon fusion contribution. However even 
so the WBF selection is composed of one-to-three parts 
GF to WBF, which means measurements of physics that 
only enters the weak boson fusion contribution will need 
to take this gluon fusion “pollution” into account. 

C. Constraining the quartic VVhh contribution 

As mentioned in Section m there is a contribution 
from quartic VVhh vertices to the WBF induced signal, 
and modifications of the corresponding gvvhh couplings 
away from their SM values using the Higgs Cross Section 














FIG. 8: Expected limits on the gauge-Higgs quartic couplings £ = gvvhh/gvvhh under the assumption of no systematic 
uncertainties (a) and 20% systematic uncertainties (b). 


Working Group k framework [57j will greatly enhance the 
signal cross section. This allows us to constrain £ defined 
by gvvhh = C x 9vvhh- To achieve this we have generated 
events with varying £ using Mad Event v5 and applied 
the WBF selections described in Section HVBI to estimate 
the enhancement of the signal, which is compared to ex¬ 
pected cross section limits on the signal with 3 ab -1 of 
data in the WBF selection under the assumptions of no 
systematic uncertainties and 20% total systematic un¬ 
certainties for comparison. The results are presented in 
Figure [8] We find that in the more realistic scenario of 
20% systematic uncertainties the expected constraint on 
the gvvhh couplings is 0.55 < C < 1-65 at 95% confi¬ 
dence level. A measurement of pp —> hhjj is therefore 
crucial to constrain new physics which enters predomi¬ 
nantly through enhancements to gvvhh- 


D. Event shapes of the tagging jets system 

The analysis strategies outlined so far have mainly re¬ 
lied on exploiting correlations in the di-Higgs system, 
with only Arj(jl.j2) carrying information about the tag¬ 
ging jets. Following similar applications in the context 
of single Higgs production |5ll| . we investigate a range 
of event shapes in the tagging jets system in the follow¬ 
ing, which could offer additional discriminating power 
through capturing colour correlations in the different sig¬ 
nal contributions beyond angular dependencies. More 
specifically, we will focus on iV-jettiness [Hi, |j63| and 
thrust major which provided the best results. 

We calculate iV-jettiness by minimising 

TN = C ^ p T ,k min( Ait^i,..., AR k , N ) (3) 

k 

where C is a normalisation which cancels when taking 


the ratio of two rs, the sum is taken over all visible mo¬ 
menta which do not belong to one of the identified Higgs 
candidates within \g\ < 5, and A Rk, n is the distance in 
the plane between the fc-th momentum and the n-th 
reference vector. r 3 / 2 is then explicitly given by T3/T2. 

Thrust major is defined by 


E k Ip* • n l 


Tjnaj = max 

n-n T =0 | Pfc | 


where is the normalised thrust vector 


117 = max 


Efc IPfc ■ n| 

E k IPfcl 


(4a) 


(4b) 


Again the sums run over all visible momenta which do not 
belong to one of the identified Higgs candidates within 

M < 5. 

We find r 3 / 2 and T ma j show promise for improving the 
WBF selection, but the signal cross section is already 
too low for us to be able to make meaningful use of this 
insight. The r 3 / 2 and T ma j distributions after the GF 
and WBF selections have been applied are presented in 
Fig. [9] Cutting, e.g., on T ma j < 0.05, the gluon fusion 
contribution is reduced by 80%, while the WBF contri¬ 
bution is reduced by only 55% amounting to a total of 2 
expected WBF and 0.3 expected GF events, with back¬ 
grounds very strongly suppressed. This means that WBF 
can in principle be observed at a small rate that can be 
used to set constraints on new physics in an almost GF- 
free selection with greatly reduced backgrounds. 

The event shape distributions can also be used to 
grea tly r educe the background in the GF selection, 
Fig. 9(c) It should be noted that these improvements of 
GF vs WBF vs background ultimately depend on under¬ 
lying event and pile up conditions and have to be taken 
with a grain of salt at this stage early in run 2. However 
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FIG. 9: Shape comparisons of A-jettiness and thrust calculated in the major direction after the gluon fusion selection of 
Sec. IIV Al (a,c) and WBF Selection of Sec. II V Bl (b,d) have been applied. 


the clear separation that can be achieved with these ob¬ 
servables indicate that an analysis employing MVA tech¬ 
niques could, at least in theory, significantly improve the 
results presented here. These techniques may also prove 
useful at a 100 TeV collider where the dihiggs production 
cross-section is substantially higher 0- 

V. SUMMARY AND CONCLUSIONS 

After discovering single Higgs production at the Large 
Hadron Collider, new analysis strategies need to be ex¬ 
plored to further constrain the presence of new physics 
beyond the Standard Model. Higgs pair production is 
pivotal in this regard as constraints from multi-Higgs 


production contain complementary information, in par¬ 
ticular with respect to the Higgs boson’s self-interaction. 
Cross sections for di-Higgs production are generically 
small at the LHC, which highlights the necessity to ex¬ 
plore other viable channels than pp —> hh to enhance sen¬ 
sitivity in a combined fit at high luminosity. To this end, 
we have investigated pp —> hhjj production in detail in 
this paper. Keeping the full top and bottom mass depen¬ 
dencies, we find sensitivity of pp —> hhjj searches at the 
LHC for production in the SM and beyond. The gluon 
fusion contribution remains important at high invariant 
di-Higgs masses where the dominant backgrounds can be 
suppressed to facilitate a reasonable signal vs background 
discrimination. Unfortunately, the gluon fusion contri¬ 
bution remains large even for selections that enhance the 
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weak boson fusion fraction of pp —> hhjj events. This 
“pollution'’ is important when such selections are em¬ 
ployed to set constraints on new physics effects that enter 
in the WBF contribution exclusively. Large new physics 
effects in the WBF contribution can still be constrained, 
which we have illustrated through an investigation of the 
constraints that can be set on deviations of the quar- 
tic VVhh couplings from their SM values with the HL- 
LHC, demonstrating that a measurement of pp —> hhjj 
will provide a powerful probe of these. Employing ob¬ 
servables which are intrinsically sensitive to the differ¬ 
ent colour correlation of WBF compared to GF, the dis¬ 
crimination between GF, WBF, and background can be 
further improved. However, the signal cross section is 
typically already too small to use such a strategy to con¬ 
strain the presence of new physics if those effects are only 


a small deviation around the SM. If new physics effects 
are sizable, such an approach will remain a well-adapted 
strategy to minimise GF towards a pure WBF selection. 
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